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Abstract. Inclusive cross-sections for gluon jet production are studied numerically in the 
' perturbative QCD pomeron model for central collisions of identical nuclei at high energies. 

Two forms for the inclusive cross-sections, with and without emission from the triple pomeron 
vertex, are compared. The difference was found to reduce to a numerical factor 0.8 -j- 0.9 for 
\^ ' momenta below the saturation momentum Q s . Above Q s no difference was found at all. The 

gluon spectrum was found to be ~ A at momenta k below Q s and ~ A above it. At large 
k the spectrum goes like l/A; 2 ' 7 " 3 3 flattening with energy. The multiplicities turned out to be 
proportional to A with a good precision. Their absolute values are high and grow rapidly with 
energy in accordance with the high value of the BFKL intercept. 
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1 Introduction 

> ■ 

In view of the new experimental data on heavy-ion collisions at RHIC and future such data to be 
obtained at LHC one would like to have predictions for the spectra of produced secondaries based 
on the fundamental theory and not purely phenomenological. At present the only candidate for 
this is the hard pomeron model derived from perturbative QCD. Originally constructed for the 
description of high-energy low -x hadronic scattering (the BFKL model [1]) it has subsequently 
been generalized to hadronic or deep inelastic scattering on nuclei [2, 3] and nucleus-nucleus 
scattering [4]. The model suffers from a serious drawback related to the use of fixed and not 
running strong coupling constant. Curing it does not look too promising, since due to absence 
of ordering of momenta in the model, it also means solving the confinement problem. However 
in spite of this defect the model seems to describe high-energy phenomena in a qualitatively 
reasonable manner. Also attempts to include the running of the coupling in some effective way 
have shown that the effect of the running is not at all overwhelming, although introduces some 
quantitative changes into the predictions. So, also for lack of something better, the perturbative 
QCD pomeron model appears to give a reasonable basis for the discussion of particle production 
in high-energy heavy-ion collisions. Of course due to the perturbative character of the model 
it can only give predictions for production of jets, leaving jet-to-hadrons conversion to non- 
perturbative fragmentation mechanism. 

From the start it has to be stressed that heavy-ion collision amplitudes are described in the 
model by complicated equations, whose solution is quite difficult to obtain even numerically 
(see [5] for partial results). Happily, as was shown in [6], due to Abramovsky-Gribov-Kancheli 
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(AGK) cancellations [7], to find the single inclusive distributions one does not have to solve 
these equations, but only to sum the appropriate sets of fan diagrams, which is accomplished 
by the non-linear evolution equation of [2, 3]. Still this operation involves a numerical study 
of considerable complexity So up to now there has been no consistent calculation of the jet 
spectra for realistic nuclei, although some preliminary attempts has been done in [8, 9, 10, 11]. 
In all cases however the authors relied on very drastic simplifications from the start choosing for 
the nuclear structure and/or for the gluon distributions in the colliding nuclei some primitive 
explicit forms in accordance with their own taste and prejudice. In fact these forms appear 
to be rather far from realistic ones, which correspond to actual participants and follow from 
the calculations. This gave us motivation to calculate numerically the jet spectra in heavy ion 
collisions as predicted by the hard pomeron model in a consistent manner. 

Another goal of the present calculations has been to compare the results obtained on the ba- 
sis of the expression for the inclusive cross-section which follows from the AGK rules applied to 
the diagrams with QCD pomerons interacting via the three-pomeron coupling [6] with a some- 
what different expression obtained from the colour dipole picture [12]. Our calculations show 
that these two formally different expressions lead to completely identical results at momenta 
of the order or higher than the value of the so-called saturation momentum Q s . At momenta 
substantially lower than Q s the colour dipole cross-sections differ from the ones from the AGK 
rules by a universal constant factor ~ 0.8 -r- 0.9. 

In both cases the spectra at momenta below Q s are found to be proportional to the number of 
participants (oc A for collisions of identical nuclei) and not to the number of collisions(oc A 4 / 3 ). 
Since Q s grows with energy very fast, the region where the spectra are oc A extends with energy 
to include all momenta of interest. At momenta greater than Q s the spectra grow faster than 
A but still much slowlier than A 4 ^ 3 (a numerical fit gives something like oc A 1 ' 1 ). 

Note that in the last years a few more phenomenologically oriented studies of particle produc- 
tion in nucleus-nucleus production have been presented, in the framework of the color-condensate 
model [13] solved in the classical approximation on the lattice [10] and in the saturation model 
[11]. In both approaches quantum evolution of the nuclear gluon density was neglected and the 
saturation momentum was introduced as a parameter fitted to the experimental data at RHIC. 
Although some of their predictions (proportionality of the multiplicity to A modulo logarithms) 
agree with our calculations with full quantum evolutions, the quantitative results are rather 
different. We postpone a more detailed discussion of this point until our Conclusions. 



2 Basic equations 

Our basic quantity will be the inclusive cross-section Iab(v, k) to produce a jet with the trans- 
verse momentum k at rapidity y in a collision of two nuclei with atomic numbers A and B: 

(27r) 2 da 

It can be represented as an integral over the impact parameter 6: 

I A B(y,k) = J d 2 M AB (y,k,b). (2) 

Our study will be restricted to the inclusive cross-sections at fixed impact parameter 6 = 
(central collisions). We shall also limit ourselves to collisions of identical nuclei A = B and for 
brevity denote Iaa = I A and so on. The corresponding multiplicity at fixed rapidity y and 6 = 
will be given by 

1 f d 2 k 
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where a A (b) is the total inelastic cross-section for the collision of two identical nuclei at fixed 
impact parameter b. For heavy nuclei one expects that a A (b = 0) ~ 1, so that the multiplicity 
is just the integral of the inclusive cross-section over the momenta. 

As argued in [4], in the perturbative QCD with a large number of colours the nucleus- 
nucleus interaction is described by a set of tree diagrams constructed with BFKL pomeron 
Green functions and triple pomeron vertices for their splitting and fusing. The structure of the 
interaction at the vertex is illustrated in Fig. 1, in which horizontal lines correspond to real 
gluons produced in the intermediate states and vertical and inclined lines describe propagating 
reggeized gluons. From this structure one sees that the produced gluons are contained in the 
intermediate states of the interacting pomerons, so that to get the inclusive cross-section one has 
to "open" these pomerons, that is to fix the momentum of one of the intermediate real gluons in 
them. A similar production mechanism in the old-fashioned local pomeron model was proven to 
lead to the inclusive cross-section given by a convolution of two sets of fan diagrams connecting 
the emitted particle to the two nuclei times the vertex for the emission (Fig. 2a). The proof 
was based on the AGK rules appropriately adjusted for the triple pomeron interaction [14]. It 
was later shown in [15] that the AGK rules are fulfilled for interacting BFKL pomerons. So the 
same arguments as in [14] allow to demonstrate that for the collision of two nuclei the inclusive 
cross-section will be given by the same Fig. 2a, that is, apart from the emission vertex, by the 
convolution of two sums of fan diagrams, constructed of BFKL pomerons and triple pomeron 
verteces, propagating from the emitted particle towards the two nuclei [6]. 

Taking into account the form of the emission vertex (see [6]) we obtain in our case (A = B 
and fixed 6 = 0) 

lA(y,k) = J d 2 f3d 2 re lkr [A<S> A (Y-y,r,f3)][^A(y,r,f3)], (4) 

Here <&(y, r,(3) is the sum of all fan diagrams connecting the pomeron at rapidity y and of the 
transverse dimension r with the colliding nuclei at distance (5 from their centers. One of the nuclei 
is assumed to be at rest and the other at the overall rapidity Y. A's are the two-dimensional 
Laplacians applied to <&'s. 

Later from the colour dipole formalism a slightly different form for the inclusive cross-section 
was derived in [12]. For the dipole-nucleus scattering case it corresponds to changing 

2<Z> A (y, 0, r) - 2$ A (y, (3, r) - & A (y, (3, r). (5) 

Note that in [12] it was erroneously stated that the change was from the "quark dipole" to 
the "gluon dipole" 2$ — <£ 2 . As seen from ( 5 ) it is not. In fact the change is equivalent to 
adding to the AGK contribution ( 4 ) a new one which has the meaning of the emission of the 
gluon from the triple pomeron vertex itself. Such a contribution is not prohibited in principle. 
From our point of view, taking into account the structure of the vertex shown in Fig. 1, its 
appearance is difficult to understand. However in this paper we do not pretend to discuss the 
validity of the two proposed formulas for the inclusive cross-sections on the fundamental level. 
Rather we shall compare the cross-sections which follow from them after numerical calculations. 

For the nucleus-nucleus case the recipe of [12] implies taking into account two new diagrams 
for the inclusive cross-sections shown in Fig. 2 b and c. As a result one finds, instead of ( 4 ), 
the Kovchegov-Tuchin (KT) cross-section 

2A<J? A (Y-y,r,/3)A<D A (y,r,/3)- 



lF(v,k) = ^ d 2 0re^ 



A<S> A (Y - y, r, (3)A<!> A (y, r, 0)-, A<S>\{Y - y, r, (3)A$ A (y, r, (3) 



(6) 
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Function 4> A (y,r,P) = <& A (y,r, P)/(2irr 2 ), in the momentum space, satisfies the well-known 
non- linear equation [2, 3] 

— = -H<f> A (y,q,P) ~ <j) A {y,q,f3), (7) 

where y = ay, a = a s N c /ir, a s and N c are the strong coupling constant and the number of 
colours, respectively, and H is the BFKL Hamiltonian. Eq. ( 7 ) has to be solved with an initial 
condition at y = determined by the colour dipole distribution in the nucleon smeared by the 
profile function of the nucleus. Both cross-sections ( 4 ) and ( 6 ) can be expressed via function 

h A (y,q,P)=q 2 V 2 q( f> A (y,q,P), (8) 

which has the meaning of internal gluon density in each of the colliding nuclei. One easily 
obtains for ( 4 ) 

I A (v,k) = ^l Jd 2 pd 2 qh A (Y-y,k-q,P)h A (y,q,P). (9) 

For ( 6 ) one also obtains a factorized expression similar to ( 9 ) 

IA T (V, k) = J d 2 pd 2 qh A (Y -y,k-q,(3) [w A (y, q, (3) - h A (y, q, (3)] , (10) 

where w A (y,q,P) is a new function, which however can be expressed via h A : 

k 2 f d 2 q 

w A (y,k,P) = — J - ^ 2 h A (y, k - q, (3)h A {y, q, f3). (11) 
Function h A (y,k,{3) has a normalization property [6] 



/ 



^h A (y,k,P) = l (12) 



and at sufficiently high y acquires a scaling property 

h A (y,k(3) = h(k/Q s (y,P)), (13) 

where Q s (y,f3) is the above-mentioned saturation momentum. From ( 12 ) and ( 13 ) one easily 
establishes some properties of the new function w A . Obviously it scales with the same saturation 
momentum when h A does 

w A (y,kf3)=w(k/Q s (y,f3)). (14) 

At k — > oo it has the asymptotic 

w A {y,k,P) k ^oo ~ 2h A (y, k,/3) (15) 

and finally 

J d 2 kw A (y, k, 0) = 2 J d 2 kh A (y, k, 0). (16) 

These properties immediately allow to make some preliminary comparison between the cross- 
sections given by ( 4 ) and ( 6 ). Obviously if k/Q s is large both expressions give the same 
cross-section due to ( 16 ). In the opposite limit of small k/Q s , the scaling property allows to 
conclude that the ratio of the two cross-sections is a universal constant which does not depend 
on y, nor on A nor on p. Our numerical results confirm these predictions. 
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3 Results 

In our study we have taken the initial condition in accordance with the Golec-Biernat distribution 
[16], duly generalized for the nucleus: 

^•'•« = 4^« Ei (-a2»)' <17) 

with 

a A = 20.8mbAT A {(3). (18) 

where T A {(3) is the standard nuclear profile function, which we have taken from the Woods-Saxon 
nuclear density. 

Evolving 4>A(y,q,(3) up to values y = 8 we found the inclusive cross-sections ( 4 )and ( 6 ) 
at center rapidity for energies corresponding to the overall rapidity Y = Y/a. with Y = 16. 
Taking a s = 0.2 this gives Y ~ 80. This value is is far beyond the present possibilities and 
was chosen only to follow the asymptotical behavior of the cross-sections at super-high energies 
and compare it with the situation at energies available presently or the near future. The overall 
cutoffs for integration momenta in Eq. ( 7 ) were taken according to 1.10" 16 GeV/c < q < 1.10+ 16 
GeV/c. For the (fixed) value of the strong coupling constant we have taken a s = 0.2 

We first discuss the cross-sections ( 4 ) obtained from the AGK rules. They are illustrated 
in Figs. 3-5. To see the absolute values of the inclusive cross-sections at different energies, in 
Fig. 3 we present them for A = 9, Y = 4, 8, 16 and y = Y/2. To illustrate the change of their 
form with energy, in Fig. 4 we present the same distributions at k > 0.3 GeV/c normalized to 
unity and multiplied by k 2 to exclude the trivial 1/k 2 dependence present in ( 4 ), 

J 9 (y,k) = k 2 I 9 (y,k)/ J ^I 9 ( y ,k'), k' > 0.3GeV/c (19) 

again at Y = 4, 8, 16. One clearly observes that below a certain point all the momentum 
dependence is reduced to the trivial factor 1/k 2 , implying that the integral in Eq. ( 4 ) is 
independent of the momentum. This point roughly coincides with the saturation momentum 
Q s (y,(3). However one should have in mind that for a given nucleus the value of the saturation 
momentum varies depending on the nuclear transverse density at distance (5 from the center 
of the nucleus. For A > 9 and values of [3 inside the nucleus we find Q s ~ 20 4- 200 GeV/c 
at Y/2 = 4 and Q s ~ 1.5. 10 5 4- 1.5. 10 6 at Y/2 = 8. Comparing with Fig. 4 we find that at 
k < Q s the integral factor in ( 4 ) is practically a constant. For k > Q s it rapidly falls. At large 
k » Q s the inclusive cross-sections at the center y = Y/2 are found to behave like l/k p (y) 
with power p(y) diminishing with energy. From our calculations we approximately find that 
p(y) = 3.3, , 3.0 and 2.7at y = 2, 4 and 8 respectively. At infinite energies p seems to tend to 2 
in correspondence to the fact that Q s — ► oo 

In Fig. 5 we illustrate the yl-dependence showing ratios 

RP art = 9lA(^ 

A AI 9 (y,k) 

with inclusive cross-sections scaled by A, at y = Y/2 and Y = 4, 8 and 16 (from top to bot- 
tom) . One clearly sees that whereas at relatively small momenta the inclusive cross-sections are 
proportional to A, that is to the number of participants, at larger momenta they grow with A 
faster, however noticeably slowlier than the number of collisions, roughly as A 1 ' 1 . The interval 
of momenta for which I a oc A can also be related to the value of the saturation momentum 
Q s (y,/3). Inspection of Fig. 5 shows that the distributions are proportional to A at values of k 
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smaller or in the vicinity of the value of the saturation momentum. Since Q s {y,(3) grows with 
energy, one may conjecture that at infinite energies all the spectrum will be proportional to A. 

Passing to the determination of multiplicities one has to observe certain care because of the 
properties of the perturbative QCD solution in the leading approximation embodied in Eqs. (4) 
and (7). As follows from ( 4 ) the inclusive cross-section blows up at k 2 — > independently 
of rapidity y. So the corresponding total multiplicity diverges logarithmically. However, the 
physical sense has only emission of jets with high enough transverse momenta. Thus one has 
to cut the spectrum from below by some krnin which separates the spectrum of jets proper 
from soft gluons which are not related to jets. Inevitably the multiplicity of thus defined jets 
depends on the chosen value of k m i n . We have chosen k m i n = 2 GeV/c. At all energies the 
multiplicities at the center ha(h = Y/2) were found to be approximately proportional to A. The 
ratios ha(v = Y/2)/A are presented in Fig. 6. One observes that their values are quite high 
and grow very fast with energies. This is not surprising, considering a very high value of the 
pomeron intercept in the lowest order BFKL model. The bulk of the multiplicity comes from 
jets with relatively high momenta. To illustrate this point in Fig. 7 we show the dependence of 
the central multiplicity on k m i n for A = 9 and at Y = 8 in the interval k m i n = 0.3 -j- 16 GeV/c. 
One observes that it goes down very slowly, indicating that it is the high momentum tail of the 
distribution which matters. 

Finally we pass to the cross-sections obtained with the KT formula ( 6 ). In Fig. 8 we show 
the ratios of these cross-sections to the ones defined by the AGK rules, Eq. ( 4 ), for y = Y/2 
and Y = 4, 8 and 16. These ratios turn to unity at k in the vicinity and above Q s , as discussed 
in the end of the preceding section. Below Q s the ratios are approximately equal to 0.8 0.9 
with little dependence on A and Y. Some dependence which is left can be explained by the fact 
that for very peripheral parts of the nucleus the scaling regime can only be reached at rapidities 
well above the considered ones. Due to this very simple relation between the two cross-sections, 
all conclusions about the A -dependence drawn for the AGK cross-section ( 4 ) remain valid also 
for the KT cross-section ( 6 ). 

4 Conclusions 

We have calculated the inclusive cross-sections for gluon production at mid-rapidity in nucleus- 
nucleus central collisions in the perturbative QCD approach with a large number of colors. 
Realistic nuclear densities were employed to account for the peripheral parts of the nuclei, whose 
contribution rapidly grows with energy due to smallness of unitarizing non-linear effects. The 
form of the cross-sections is found to be determined by the value of the saturation momentum 
Q s , which depends on the rapidity and nuclear density. At momenta much lower than Q s 
the spectrum is proportional to 1/k 2 . Its A dependence is close to linear. At momenta much 
higher than Q s the spectrum is found to fall approximately as l//c 2 - 7 ~ 3 - 3 with the yl-dependence 
as ~ A 11 . The multiplicities at mid-rapidity are found to be proportional to A with a good 
precision. They grow with energy very fast which is related to a fast growth of the saturation 
momentum. 

We also compared two different forms for the inclusive cross-section, which follow from the 
AGK rules or the dipole picture. The difference between their predictions was found to be absent 
for values of momenta larger than Q s . At momenta smaller than Q s the difference reduces to a 
universal constant factor: the dipole cross-sections are just 0.8 -j- 0.9 of the AGK cross-sections. 
With the growth of Y and/or A this factor slowly grows towards unity, so that it is not excluded 
that at infinite A and Y the two cross-sections ( 4 ) and ( 6 ) totally coincide at all values of 
momenta. All our conclusions about the energy, momentum and A dependence are equally valid 
for both forms of the inclusive cross-sections. 

As mentioned in the Introduction a few more phenomenological studies of the gluon pro- 
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duction in nucleus-nucleus collisions were recently made in the classical approximation to the 
colour-glass condensate model [10] and in the saturation model of [11]. In both studies quan- 
tum evolution was neglected, so that scaling with the saturation momentum Q s was postulated 
rather than derived. The saturation momentum thus appeared as an external parameter, whose 
A and Y dependence were chosen on general grounds and whose values were fitted to the exper- 
imental data at RHIC. In both models the multiplicities turned out to be proportional to the 
number of participants (modulo logarithmic dependence on A different in the two approaches). 
This agrees with our results. However the form of the inclusive distributions in momenta found 
in [10] is different from ours. Its behaviour both at small k (~ 1/y/k 2 + m 2 with m = 0.0358Q S ) 
and at large k (~ l/k A ) disagrees with the form of the spectrum we have found. For realistic 
nuclei the spectrum was calculated in [10] only up to 6-^7 GeV/c, so it is not possible to see 
if any change in its A-behaviour will occur at higher momenta. But most of all, the value of 
the saturation momentum and the speed of its growth with rapidity which we have found from 
the QCD pomeron model with full quantum evolution are much larger than the fitted values in 
both [10] and [11]. This is no wonder in view of a very high value of the BFKL intercept in 
the leading approximation. From the phenomenological point of view this is the main drawback 
of the BFKL theory. To cure it one possibly has to include higher orders of the perturbation 
expansion and the running coupling constant. Although some work in this direction has been 
done for linear evolution [17], no attempts to generalize this to non-linear evolution in some 
rigour has been made yet. 
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Figure 1: Interaction of three BFKL pomerons at the spliiting vertex 
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Figure 2: Typical diagrams for the inclusive cross-section in nucleus-nucleus collisions. 
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Figure 3: Inclusive cross-sections Ig(y,k) at y = ^ . Curves from top to bottom at small k 
correspond to scaled overall rapidities Y = 4, 8, 16. 
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Figure 4: Normalized distributions Jg(y,k) (Eq. (25)) at y = ^ . Curves from top to bottom 
at small k correspond to scaled overall rapidities Y = 4, 8, 16. 
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Figure 5: From top to bottom: A-dependence of momentum distributions, scaled with the num- 
ber of participants, at Y = 4, 8, 16. Curves from bottom to top show ratios 9Ia(u, k)/AIg(y, k) 
at center rapidity (y = \) for A = 9, 27, 64, 108 and 180. 
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Figure 6: From top to bottom: A-dependence of multiplicities, scaled with the number of 
participants, at Y = 4,8,16. Curves from bottom to top show fiA(y)/A at center rapidity 
( y = Y) for A = 9, 27, 64, 108 and 180. 
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The dependence of multiplicity for A = 9 at center rapidity and Y = 8 on the infrared 
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Figure 8: Ratios of the KT inclusive cross-sections, Eq. ( 6 ), to the ones found from the AGK 
rules, Eq. (4) at center rapidity (y = Y/2) and Y = 4,8,16. Curves from bottom to top refer 
to A = 9, 27, 64, 108 and 180 



